Random wave functions and percolation 



E Bogomolny and C Schmit 



CNRS, Universite Paris-Sud, 

Laboratoire de Physique Theorique et Modeles Statistiques, 
91405 Orsay Cedex, France 

E-mail: eugene .bogomolnyOlptms .u-psud.f r 

Abstract. Recently it was conjectured that nodal domains of random wave functions 
are adequately described by critical percolation theory. In this paper we strengthen 
this conjecture in two respects. First, we show that, though wave function correlations 
decay slowly, a careful use of Harris' criterion confirms that these correlations 
are unessential and nodal domains of random wave functions belong to the same 
universality class as non-correlated critical percolation. Second, we argue that level 
domains of random wave functions are described by the non-critical percolation model. 
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In 1977 Berry pj conjectured that wave functions of generic chaotic systems, can 
statistically be written as random superposition of a complete set of functions, ipm{x), 

*(f) = 5]C„V^^(f) (1) 

m 

where the coefficients Cm are independent identically distributed random variables with 
zero mean and variance obtained from normalization. In particular, any wave function 
of a two-dimensional billiard obeys the Helmholtz equation with energy E = k"^ 

(A + E)*(x,i/) = (2) 

and can be represented as a formal sum 

^{x,y) = Y.CmJ\m\{kry'^^ , (3) 

m 

where r and are polar coordinates and Jm{f) are the usual Bessel functions. For 
problems without magnetic field wave functions are real and Cm = C*_^. 

Consider a chaotic quantum system like the stadium billiard with Dirichlet 
boundary conditions and let us calculate (numerically) a large number of eigenfunctions 
with energies close to fc^. Each eigenfunction gives a well defined set of coefficients. Cm 
in the expansion ([3]) and it is of interest to know mean values of different functions of 
these coefficients over the whole ensemble of eigenfunctions. Berry's conjecture means 
that in the semiclassical limit k ^ oo the result will be the same as if the coefficients 
Cm would be independent Gaussian random variables with 

{Cm)=0 and {CmC:) = a^6mn . (4) 

Though Berry's conjecture is one of the oldest conjectures in quantum chaos it has not 
yet been proved rigorously. Numerically it works very well. In figure [1] we present two 
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pictures, the left one is a true eigenfunction of a quarter of the stadium with area equal 
An and the right is a random realization of ([3]) with approximately the same energy. 
In these figures back and white regions represent nodal domains where the function is, 
respectively, positive and negative. The two figures look very similar and the left figure 
was magnified in order to distinguish better the circular arc of the stadium billiard 
where the true wave function is zero. 

In 2002 Blum, Gnutzmann and Smilansky attracted wide attention to such type of 
pictures by initiating the investigation of nodal domains of wave functions for different 
systems [2]. In [3] it was argued that nodal domains of chaotic wave functions (like in 
figure [I]) can adequately be described by a critical percolation model. This conjecture 
gives a very detailed description of nodal domains and a large number of interesting 
(global) quantities can directly be transposed from percolation results (see e.g. [4]). 

In particular, it was shown in [3j that the total number of connected nodal domains 
for random wave functions has Gaussian distribution with mean value, n{E), and 
variance, a'^{E), given by the following expressions 

n{E) 3^3-5 



and 
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Here N{E) is the mean number of eigenvalues with energy less than E. For billiards, 
N{E) ^ Ak'^/An where A is the billiard area. 

In [3] it was demonstrated that these relations are well fulfilled by random wave 
functions. For completeness we plot in figure [2] the energy dependence of number of 
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Figure 2. Number of nodal clusters for stadium eigenfunctions vs energy. Solid line 
is the asymptotic prediction ([5]). 
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Figure 3. Left: The stadium eigenfunction with E = 10098.531. Right: The same 
but with E ^ 10107.147. 



nodal domains computed numerically for true eigenf unctions of a quarter of the stadium 
close to the 10000*^ level. This fi gure confirms that expression ([5]) agrees well with the 
results of numerical calculations for chaotic wave functions. Still the numerical data 
present deviations from the expected value ([5]) larger than it follows from ([6]). The 
origin of such large fluctuations is related with the existence of wave function scars. For 
illustration we present in figure [3] pictures of nodal domains corresponding to the two 
lowest points in figure O It is clearly seen that these pictures have two different parts. 
One shows a regular pattern typical for integrable systems and the other is practically 
indistinguishable from a random wave function as in figure [H It is evident that quasi- 
regular regions are related to the bouncing ball scar. The left eigenfunction corresponds 
to the single excitation of the bouncing ball and the right one demonstrates the double 
perpendicular excitation. Regions of quasi-integrable behavior have considerably less 
number of nodal domains than chaotic regions which explains the wide distribution of 
points in figure O 

In [3] it was also checked that the distribution of nodal domain areas agrees with 
the percolation theory prediction 

n(s) ~ s'^ with T = 187/91 (7) 

and the the fractal dimension of nodal clusters is close to the percolation value 

D = 91/48 . (8) 

These and other numerical calculations confirm that nodal domains of random 
wave functions (and of chaotic quantum systems) are correctly described by critical 
percolation theory. 

However, this conjecture has an intrinsic difficulty [5]. The standard percolation 
theory deals with the following situation. One has a lattice with each site (or edge) 
positive with probability p or negative with probability 1—p. The important assumption 
here is that these probabilities (or concentrations) at different points are independent 
random variables. For random and chaotic wave functions this is not the case. From 
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([2]) it follows that 

{^{x,y)) = (9) 
and the values of the wave function at two points are correlated 

{<i!{xi)^{x2)) = a^g{\x,-X2\) (10) 

where 

oo 

^7(15^1 -X2\)= E Jm{kri)Jmikr2)e'"'^'''-'''^ = Mk\x^ - fsl) • (H) 

m=~oo 

In mathematical literature one often defines a Gaussian random function without explicit 
representation as in ([1]) but by postulating its mean value and variance (cf. Qj and ffTOj) ). 
all other correlation functions being determined by the Wick theorem. 

It is well known [6] that for a Gaussian random function the probability that its 
value at two points has the same sign is 

P(|fi -fsl) = I + -arcsinG(|xi -fal) (12) 

Z TT 

where G{\xi — X2I) is the two-point correlation function divided by its value at so that 
G{0) = 1. Therefore, if a random Gaussian function is positive at a point, at all close-by 
points it will be also positive with probability close to 1. Only at large distances when 
G{r) — > probabihties of being positive and negative become equal as in the usual 
percolation. For random wave functions G{r) coincides with g{r) in ( iTTl) and because 



Mkr) ™ \^cos( kr - ttM) (13) 
V nkr 

the random wave function correlation function decays quite slowly (~ r~^/^) which may 
destroy the validity of percolation theory. 

One of the purpose of this paper is to show that this is not the case. In section [2] 
we use so-called Harris' criterion [7j to demonstrate that the oscillating character of 
the correlation function ffTTl) leads to strong cancelations and even the slow decay as 
in (|T3l) is sufficient in order that nodal domains of random wave functions be in the 
same universality class as non-correlated percolation. Though Harris' criterion is not 
mathematically rigorous, it is widely accepted (see e.g. [8j) and its fulfillment is a 
strong argument in favor of the validity of the conjecture that nodal domains of random 
wave functions in the semiclassical limit k 00 are described by critical percolation. 
In section [3] we generalize this conjecture to level domains i.e. regions where \l/(a:, ?/) 
is bigger (or smaller) than a non-zero value e. We argue that such domains are also 
described by percolation theory but this time by the non-critical percolation and the 
deviation from criticality is proportional to e. 



2. Harris' criterion 



Harris proposed his criterion in 1974 [7] to investigate the infiuence of random defects on 
the critical behavior of the Ising model. Later it was applied to many other models. In 
particular, in ^ it was used for correlated percolation and we follow closely this paper. 
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Let us consider a critical percolation problem where concentrations in different 
sites, denoted by ^(x), are correlated. The criticality condition means that the mean 
concentration in each site equals a critical value p* 

{^{x))=p*. (14) 

Below we assume that the connected part of the two-point correlation function of 
concentrations depends only on the difference between these points 

((^(fi) ~p*){'^{x2)-p*)) = 9{\xi - X2I) . (15) 

Harris' criterion was developed to give a condition on function g{r) under which 
correlations are unessential and critical properties of the problem under consideration 
are the same as for non-correlated percolation. Harris' argumentation is as follows [3 [8]. 
Let pv be the average concentration in a finite volume V 

pv = hi: *(^) • (16) 

The mean value of pv over realizations equals the critical value 

{pv)=P* (17) 

but for a given realization it differs from it. 

Consider the situation when inside the volume V all sites have the same 
concentration pv- As pv 7^ p*, it corresponds to a non-critical percolation. From 
the percolation theory it is known that points are uncorrelated if they are separated by 
a distance larger than a certain correlation length (the cluster size) ^ 

where u is certain critical index. For the standard two-dimensional percolation z/ = 4/3 
When py is close to p* this correlation length is large but finite. Divide the whole 
space into cells of radius of order of ^ in (fTSl) . Different cells can be considered as 
uncorrelated but inside each cell concentrations are correlated. The natural measure of 
the importance of correlations is the variance, A, of pv which can be computed from 
the knowledge of the two-point correlation function (fT5|) 



{{py-pr) = j-, E m^i) - p*mx2) - p*)) 



g{\xi - f2|)dxidf2 • (19) 

xiGV Jx2&V 



1 

Harris' criterion [7], [8] states that if the variance is small with respect to \py — 

A«by-p*P, (20) 

correlations are unessential and all critical quantities are the same as for the standard 
uncorrelated percolation. 

As the integrant in ( |T9l) depends only on the distance between two points one 
usually simplifies this integral without discussion (cf. [8]) by assuming that one point 
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is somewhere inside V and the second point is on a distance r of the order of the initial 
cell size so that the double integral in f|T9l) can be rewritten as a simple one 

^^hL amdx. (21) 

When the d-dimensional cell V is of the radius ^, Eq. (pT!) leads to the following widely 
used expression 



If the correlation function decreases as a certain power of the distance 



9[r) 



one gets 



Therefore 



dr 



const when a > d 



when a < d 



(22) 
(23) 
(24) 



A 



v-P 



* I du 



(25) 



when a > d 
~ \Pv - P*]"" when a < d ' 
The comparison with ( l20l) leads to the standard Harris criterion ^ that correlations are 
unessential provided min(a,(i)z/ > 2. 

For random wave functions g{r) decays as in ( fT3l) and because for the two- 
dimensional percolation u = A/3 Harris' criterion seems to indicate that random waves 
cannot be described by usual percolation model. But the above considerations are valid 
only when the correlation function g{r) is positive. When it has oscillations, formulae 
(pTj) and (122!) cannot be correct as by definition A is non-negative but these expressions 
may change the sign. 

The main point is that the double integral in f|T9|) cannot, in general, be 
approximated by the simple integral ( l2Ti) . 

For random wave functions correlation function is given by (ITT!) 



g{\xi - X2\) ~ Mk\xi - X2I) = J\m\{kri)J\m\{kr2)e 



im(0i-(/)2) 



(26) 



m=0 



Assuming for simplicity that region y is a circle, all integrals over polar angles vanish 
except the one with m = and 

g{\xi — X2\)dxidx2 = (27r ( jQ{kr)rdr\ . 



Therefore in this case 
Finally as z/ = 4/3 one concludes that 



A 



Jo(kr)rdr 



—3 



(27) 



(2^ 



(29) 



which means that Harris' criterion is fulfilled and nodal domains of random (and chaotic) 
wave functions in the universal (continuum limit) have the same critical behavior as the 
standard critical percolation model. 
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Figure 4. Left: Nodal domains of a random wave function. Right: Level domains of 
the same fmiction with e = 0.03. In the both figures the largest connected clusters are 
highlighted. 

3. Level domains and non-critical percolation 

The success of percolation-like description of nodal domains of random wave functions 
naturally leads to different generalizations. In particular, it was mentioned in [3] that 
level domains of random functions can also be described by percolation theory. But 
contrary to the above case one has to consider not critical but non-critical percolation 
theory. 

Level domains are regions where a function is bigger than a certain value 

^{x,y)>e (30) 

with e 7^ 0. To understand better how level domains look like we present in figure H] 
two pictures. The left one displays nodal domains for one realization of a random wave 
function ([3]). Boundaries between black and white nodal regions correspond to nodal 
lines where y) = 0. The right figure represents level domains of the same wave 
function. White (resp. black) regions are domains where \l/(x,|/) is bigger (resp. lower) 
than 0.03 (in a certain normalization). Boundaries between these regions (called level 
lines) are solutions of equation \l/(a;, y) = 0.03. Both figures look quite similar. To stress 
their difference the largest connected clusters in both figures are highlighted. In the left 
figure this cluster connects different boundaries but does not cover the whole region. 
On the contrary, in the right figure the cluster fills practically all allowed space. This 
behavior is quite reminiscent of percolation transition. At critical concentration there 
exist an infinite cluster and a large number of finite clusters (cf. (I7l) and (IE])). When the 
concentration is bigger than the critical, one infinite cluster with the dimension equal 
to the space dimension appears. 
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To quantify different relations between level domains and non-critical percolation 
we use one-parameter scaling which is assumed to be valid for non-critical percolation 
not too far from criticality. According to the scaling conjecture all quantity for non- 
critical percolation with concentration p are the same as for critical percolation with 
concentration p* multiplied by a function depending only on one universal combination 
of parameters. 

For example, the number of cluster of size s in critical case is given by ([7]). For 
non-critical percolation it takes the following form 

ns{p) = K,s^f{K2{p-p*)sn (31) 

where f{z) is a certain function depending on the universal argument 

z={p-p*)s''. (32) 

Here r and a are universal critical indices. For two-dimensional percolation they are 
known analytically (see e.g. |4j) 
187 36 

T = , a = — . 33 

91 91 ^ ^ 

In fl3T]) only the constants Ki, K2, and the threshold value p* depend on micro details 

of the problem such as the type of percolation (side or edge percolation) and the form of 

the lattice (triangular, square, etc.). The indices r, a and the function f{z) are assumed 

to be universal (but depend on dimensionality of the problem). 

Scaling and universality are very important properties of non-critical percolation 

theory. We will now check that they are also well fulfilled for level domains of random 

wave functions. First we need a relation between the deviation from criticality, p — p*, 

and e in (!30|) . To simplify the formulas below it is convenient to normalize billiard wave 

functions by the condition 

\'^{x,y)\'^dxdy = A (34) 

A 

where A is the billiard area. 

With this normalization random wave functions are Gaussian random functions 
with (7 = 1 (cf. (fTTj) ). Therefore, the probability density that "^{x,!/) = t is 

Pit) = (35) 
V 27r 

and the probability P>{s) that \E'(x, y) > e takes the form 

P>{e) = I P{t)dt ^ 2 - + • (36) 

Consequently, for small e 

p-Pc-^ e (37) 

and for level domains the scaling relations as in (13T!) are valid with the substitution 

p — p* — > e. 
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Figure 5. Ratio of the number of clusters to its critical value for level domains of 
random wave functions (solid line) and for non-critical bond percolation on a square 
lattice. For the first curve the abscissa axis is e, for the second one it represents p — p* . 
Black circles are results of rescaling p — p* = e/1.62. 



Consider first the total number of connected level domains. For this quantity there 
exists no explicit scale as in (13T|) and the scaling prediction is particular simple. For the 
non-critical percolation 

^-np-p') (38) 

so for level domains it should be 

W) ^ '''' 

with a certain constant K and the same function f{z) . Here n{p) is the mean number 
of connected clusters for a percolation model with concentration p. n{e) is the same 
quantity but for level domains (J30il . 77,(0) is the mean number of nodal domains given 
by©. 

In figure [5] there are two curves. The solid line is the ratio ( 139|) computed 
numerically for random Gaussian functions. Data are well fitted by a quadratic fit 
1 -|- .61356: -|- 13.47e^. The dashed line represents the ratio (138|) computed numerically 
for bond percolation on 110 x 110 square lattice. Black circles result from the rescaling 
p — p* = e/1.62. It is clearly seen that points for level domains after rescaling agree very 
well with data from pure percolation which confirms the scaling ( 138|) and ( l39l) . 
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ln(s) 



Figure 6. Number of domains of area s for random wave functions. Different curves 
represent results for different values of e. Straight line has the slope —187/91 predicted 
for nodal domains. 

In figures [6] -[7] we present the verification of fl3Tl) for level domains of random wave 
functions. In figure [6] the number of clusters with area s is plotted for a few different 
values of e and k = 200. Curves are clearly different and slopes for e ^ differ from 
the critical percolation value predicted for nodal domains. Then we rescale these 
data using as the abscissa axis the variable z = es'^ as it follows from the one-parameter 
scaling. The resulting points are concentrated close to each other and in the left of 
figure [7| the solid line indicates the average of all such calculations which supports the 
existence of scaling for level domains. In the same figure two more curves are presented. 
The dashed line represents the results of similar calculations but for the non-critical 
bond percolation on a square 110 x 110 lattice and the dotted line is the same but on 
a square lattice with 50 x 50 sites. For these curves the z variable is defined as in ( pTl) . 
namely, z = {p — p*)s'^. Once more the scaling relation is well satisfied. Finally on the 
right part of figure [8] it is shown that after the rescaling p — p* = e/1.33, the curve of 
level domains practically coincides with the one for non-critical percolation on a square 
50 X lattice. 

Finally in figure M the results of numerical calculations for fractal dimensions of 
clusters are presented. According to percolation theory (see e.g. [4J) the dominant 
cluster has the dimension 2 but the other ones have the fractal dimension equal 
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Figure 7. Left: Scaling relations. Solid line - random wave functions with k = 200. 

Dashed line bond percolation on a square 110 x 110 lattice. Dotted line - the same 
but on a 50 X 50 lattice. Right: Rescaling of dotted curve hy p — p* = e/1.33. 




Figure 8. Fractal dimension of level domains. From bottom to top: e = 0, e = 0.1, 
and e = 0.2. For all curves k = 200. The slope of the solid line is 91/48, of the dotted 
line is 2, of the dashed hue it is 1.56. Black circles represent domains where \l/(a;, y) < e 
and open ones regions with ^{x,y) > e. 



approximately to 1.56. Notice that for critical percolation (and the nodal domains) 
the fractal dimension of clusters is 91/48. 
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4. Conclusion 

The conjecture that random (and chaotic) wave functions in the semiclassical limit can 
be described by percolation theory is analyzed. We demonstrate that, though wave 
correlations decay slowly, a careful application of Harris' criterion suggests that nodal 
domains of random wave functions are indeed in the same universality class as the 
critical non-correlated percolation. We also present arguments in favor that the non- 
critical percolation is useful in the investigation of non-zero level domains of random 
wave functions. Let us mention that in [9] it is shown that nodal lines of random wave 
functions are close to SLEg curves which describe boundaries of critical percolation 
clusters. All these results confirm the above conjecture. It is of interest to find a 
rigorous proof of it. 
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